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Abstract 

A typical problem of the leptogenesis scenario is the mismatch between 
the maximum reheat temperature implied by gravitino overproduction bound 
and the minimum temperature required to create thermally the lightest right- 
handed neutrino. We explore the possibility of baryogenesis via leptogenesis in 
the presence of low scale mass right-handed neutrino. In such a scenario, right- 
handed neutrinos are created thermally at low reheat temperatures without 
relying on non-perturbative production mechanisms. We focus on two specific 
realizations of the scenario, namely the out-of-equilibrium decay of right-handed 
neutrinos (Fukugita-Yanagida) and the leptogenesis via the LH U flat direction 
(Affleck-Dine). We find that in general, the two scenarios are able to produce 
the required baryon excess for a reasonable amount of CP violation. 
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1 Introduction 



Recent experimental results gave overwhelming evidence that neutrinos have small 
but non- vanishing masses [ffl. In the standard model (SM), neutrinos are exactly 
massless, hence the explanation of neutrino experiments requires Physics beyond the 
SM. Furthermore, neutrino masses appear to be very small with respect to the other 
fermions ones. If neutrinos are Majorana particles, it is possible to accommodate 
small neutrinos masses in the SM by introducing the lepton number violating effective 
operator O e s = aij£jT2f£jH T T2fH/M, where ^ and H are the lepton and the 
Higgs doublet respectively. Here M is the scale where "new Physics" is expected 
to occur, is usually taken as the Planck or the GUT scale. In the former case, the 
presence of this lepton number violating operator is motivated by the common belief 
that gravity does not respect any global quantum number 0, |j, or at least this is 
what happens for example in black holes and wormholes -no hair theorems. In the 
latter case (M = Mqut), the effective operator arises via the see-saw mechanism 0, 
when integrating-out the heavy right-handed neutrinos (RHNs hereafter). 

On the other hand, our Universe appears to be constituted exclusively of baryons. 
In order not to spoil the Big Bang Nucleosynthesis (BBN) successful predictions 
of the observed light elements abundances ||, a small baryon excess have to be 
present. The required value is quantified by the baryon-to-entropy ratio and is given 
by Yb = tib/s = (7.2 ± 0.4) x 1CT 11 0. To accomplish successfully their task, 
baryogenesis scenarios M have to satisfy three essential conditions 0, namely: (i.) 
Baryon number violation, (ii.) C and CP violation, and (Hi.) Departure from thermal 
equilibrium. One particularly appealing scenario is the leptogenesis scenario, where 



lepton number, produced either by the out-of-equilibrium decay of heavy RHN's ||T0| 



or by the decay of a scalar condensate carrying non-zero lepton number |TT|, [12 
is reprocessed to a baryon asymmetry via the sphalerons interactions. Given the 
experimental evidence that lepton number is violated in neutrino oscillation and the 
fact that proton decay have not been observed yet, the present experimental situation 
seems to favor this scenario over the other existing baryogenesis scenarios. 

A generic problem of thermal leptogenesis scenarios is the mismatch between the 
maximum reheat temperature implied by gravitino overproduction and the minimum 
temperature required to thermally create heavy RHNs Trh 10 10 GeV. To reconcile 
these two facts, non-thermal creation of RHNs in a low reheat temperature plasma 
were considered. These mechanisms, however, involve non-perturbative dynamics 
and are in general sensitive to inflation models. Furthermore, they lead to even more 
stringent bounds on the reheat temperature, due to the non thermal production of 



moduli and gravitinos [13, 14 



The aim of this paper is to address this issue in a different perspective. We will 
consider the situation where the reheat temperature is low (may be as low as the TeV) 
and we will only consider thermal production of RHNs. This will naturally lead us 
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to consider a class of see-saw models (that we will subsequently call low-scale see-saw 
models), where RHNs have TeV masses instead of the conventional unification scale. 

The paper is organized as follows. In section |2], we give our main motivation for 
the scenario. In section p|, we study leptogenesis through the out-of-equilibrium decay 
of low scale RHNs. In section |], we turn to the Affleck-Dine scenario. Finally, in 
section || we summarize our conclusions. 



2 The gravitino problem vs. thermal leptogenesis 

As any unwanted relic, gravitinos represents a potential danger for the thermal history 
of the Universe. Gravitinos are created predominantly via 2^2 inelastic scatterings 
of gluons and gluinos quantas. Their relic density and contribution to the energy 



density are given by |T5 



v 1 i x in-" ( Trh \ f 100 GeV V ( m ' 9 \ 



J 



where rrig denotes the gluino mass. The requirement that, if unstable, their late 
decay do not disrupt the successful BBN predictions, and if stable, their energy 
density do not overdose the Universe, put tight constraints on their relic abundance. 
It has been noted that if m 3 / 2 >10 TeV or m 3 / 2 < keV, then there is no gravitino 



problem [16], [17]]. These requirements can be relaxed if there is a period of inflation 
and the constraints apply only on post-inflation abundances. From the expression (TI]), 
one sees that the gravitino abundance scales linearly with the reheat temperature, 
therefore the bound on Y 3 / 2 translates onto the following bound on the maximum 



allowed reheat temperature Trh 18 



Trh < (10 6 - 10 9 ) GeV for m 3/2 = 100 GeV - 1 TeV . (3) 

There exists however more stringent bounds on Trh from non-thermal production. 
For generic supersymmetric inflation models, the bound can be as tight as |14[] Trh < 



10 5 (V 1//4 /10 15 GeV), where V 1 ^ is the height of the inflationary potential. 

Let us now see the constraints on the reheat temperature coming from lepto- 
genesis. In the original see-saw model || the mass scale of RHNs is typically of 
O(10 10 — 10 15 ) GeV. In addition, the bound on the CP parameter [jHJ for hierarchi- 
cal RHN's in thermal leptogenesis implies a lower bound on the mass of the lightest 
RHN M Nl > 10 10 GeV. Consequently, if the thermal lepto genesis scenario is truly 
the mechanism responsible for the the generation of the Baryon Asymmetry of the 
Universe (BAU), RHNs of this mass have to be produced after inflation. This means 
a high reheat temperature, at least as high as the mass of the lightest RHN, i.e. 
10 10 GeV, potentially conflicting with the gravitino bound discussed above. 
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A possible way out to get around this problem is to produce RHNs non-thermally, 
that is during an efficient preheating phase [p0[] . Non-thermal production, however, 



can lead in some cases to even more stringent bounds on the reheat temperature. In- 
deed, for typical hybrid inflation models, the upper bound on the reheat temperature 



can be as low as 1 TeV [13 



From the above discussion, it is clear that any compelling solution to this problem 
will, in one way or another, involve low reheat temperatures. After all, we dont know 
the thermal history of our Universe before BBN. All we know experimentally is that 
Trh > Tbbn ~ MeV. In this paper, we will consider a rather exotic solution to this 
problem, namely the case for leptogenesis when RHNs have a low scale mass. The 
first benefit of such an approach is that RHNs can be produced thermally with a low 
reheat temperature T RH ~ 0(TeV), avoiding thus the creation of dangerous relics, like 
heavy GUT monopoles, and more importantly suppressing the creation of gravitinos. 
On theoretical grounds, nothing forbids the mass of RHNs to be of 0(TeV). In fact 
this situation is encountered in many cases (See for e.g. [^1], |2^, |23], |4[]). This is also 



a typical situation that arises in models where the fundamental scale (the GUT scale 
and/or the quantum gravity scale) is of C(TeV). In this case, the Yukawa couplings 
of RHNs have to be much smaller to produce phenomenologically acceptable light 
neutrino masses. Such a fine-tuning is stable under radiative corrections because that 
Yukawa couplings are self renormalizable and is protected by supersymmetry. There 
remains the question of how such suppressed Yukawa couplings can arise in a concrete 
model. This can be achieved for example by the mean of some i?-symmetry that 
forbids the bare Yukawa coupling between the left and the right-handed neutrinos. 
As a result the leading Yukawa couplings will be suppressed by powers of a heavy 
scale [^1], |22|]. The Yukawa suppression can be obtained upon integrating-out some 



heavy field as well [25 . 



3 Thermal Leptogenesis with TeV scale RHNs 

We begin by reviewing the basics of the out-of-equilibrium decay leptogenesis scenario. 
Consider the Minimal Supersymmetric Standard Model extended by three RHNs, 
one for each generation. The interactions of the RHNs are given by the following 
superpotential 

W N = Y ij L i H u N j + hd^l (4) 

After integrating-out the RHN and electroweak symmetry breaking, the light neutri- 
nos mass matrix will be given by the familiar see-saw formula 

m u = -Y T M~ l Y (H u ) 2 . (5) 
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In this scenario, the RHNs must decay out-of-equilibrium. A measure of the departure 
from thermal equilibrium is given by the parameter K defined as 



K 



y n 

2H 



T=M N 



(6) 



where is the decay rate of RHNs and H is the expansion rate of the Universe. The 
decay is out-of-equilibrium when K < 1. The final baryon asymmetry reprocessed by 
sphalerons is given by |26| 



Y B = 



n B 



8n a + An 



22n g + 13tih J s 



(7) 



where n g and % counts the number of fermion generations and Higgses respectively. 

The lepton asymmetry produced by the CP-violating out-of-equilibrium decay of 
the RHNs can be computed using 



nL 
s 



where g* is the effective degrees of freedom and n is the dilution factor, computed by 
integrating the relevant set of Boltzmann equations |27[ [28 . The parameter e char- 
acterizing CP violation in the RHNs decay, can be defined for each RHN separately 
as @ 



Si = 



Ylj ^(Ni — > £jh u ) 



-— ^A^yim \{(yy%} 2 ] 



K\ (Mi 
M? S \M? 



(9) 



where Fy and F$ are the contributions of the vertex and self-energy respectively. 
They are given by 



Fy(x) = y/xlli 



i + i 

X 



Fs(x) 



2y/x_ 
x — 1 



(10) 



Now, applying the above formulae to TeV mass RHNs, one immediately sees that, 
due to the smallness of the Yukawa couplings, the decay of RHNs is automatically 
out-of-equilibrium. In addition to the decay processes, there can be other competing 
processes that might bring the RHNs to thermal equilibrium, depleting any pre- 
existing lepton number. These processes have to be out-of-equilibrium too, i.e. V ~ 
(nav) H. The first such process is the AL = 2 scattering £h u <-> £h u , via both 
s and t channel. Other competing processes may involve the t-(s) quark, such as 
Nt(b) <-> tb(t). It turns out that due to the Yukawa coupling suppression all these 
processes are out-of-equilibrium. Finally, it has been noted [BUI that the process 
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W /± W /± — > £ ± £ ± , mediated by virtual left-handed neutrinos can lead to stringent 
constraints on their masses. In our case, it leads to a very mild constraint. 

So far for the out-of-equilibrium conditions, now we concentrate on the CP vio- 
lation parameter e. As we have seen previously, due to the smallness of the Yukawa 
couplings, it is very easy to satisfy the out-of-equilibrium condition, however the re- 
sulting CP violation parameter e is too small. This is due to the fact that the decay 
rate V and the CP-parameter e are both proportional to the same Yukawa couplings 
combination. From Eqts (|| [10]), one sees that the two contributions to the CP pa- 
rameter e are sensible to two completely different patters of RHNs masses. While the 
vertex contribution Fy is enhanced for large hierarchies, the self-energy contribution 
Fg is so when RHNs are (quasi-)degenerate Q. In order to enhance the value of e, 
one have to exploit the properties of the two functions Fy and F$. In the next sub- 



section, we will consider the case where RHNs are nearly degenerate [32, 31]. There 



exist however another possibility, related to the fact that RHNs masses and soft SUSY 
breaking A-terms are of the same order i.e. 0(TeV). 



3.1 Leptogenesis with quasi- degenerate TeV scale RHNs 

Consider a model where two out of the three RHNs are quasi degenerate, that is iV\, 
N 2 and N 3 have masses M 1? M 2 ~ (9 (TeV) <C M 3 respectively. The mass splitting 
8M\2 = |M 2 — Mi| = 5 ■ M , where M ~ TeV). Due to their suppressed Yukawa's, 
RHNs will be long-lived enough to eventually dominate the Universe before decay- 
ing. The condition for RHNs dominance can be written Tn <C r v , where is the 



decay rate of the inflaton |j33[ . While, RHNs can hardly dominate the energy den- 
sity of the Universe because of Pauli blocking, this can happen more easily for their 
scalar partners the RH sneutrinos. Moreover, due to quantum de Sitter quantum 



fluctuations and for H in { 3> Mjv ~ O(TeV), RH sneutrinos become coherent over 
super-horizon scales and can be considered as classical fields with the constant value 
(vev) (N 2 ) = SH^/S^M^. Therefore if the RH sneutrinos scalar potential is just 
given by the mass term, they are likely to dominate quickly the energy density of 
the Universe. Given the above discussion, one can compute the lepton asymmetry 
produced during the decay of N 2 using 



til _ 3 T/y 2 
s ~ 4 Mo 



e 2 , (11) 



where T/v 2 is the decay temperature of iVVs computed by equating the energy density 
of RHNs with the energy density of the Universe when H ~ T 2 . Since N\ and iVo, are 
quasi-degenerate, we can safely ignore the vertex contribution to the CP parameter 



1 For perturbation theory to hold, the mass splitting SMik = \M{ — Mk\ must satisfy SMik ^S> T, 
where T is the decay rate of RHNs, otherwise one can no more trust the perturbative calculation 
based on Eqts (|9|,[To|) and one have to rely on a resummation approach |n| . In the limit of exact 
degeneracy, the CP parameter vanishes. 
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(Fg 3> F v ). Using Eqts (||) and (|T0D , we can compute the total CP parameter 
e ~ £i + e 2 , giving 

,s sE(pk hWo,1 l (12) 

A rough estimate of the required degeneracy gives 5 ~ 0(1CT 6 — 1CT 7 ), and pertur- 
bativity is clearly satisfied (See footnote [1] on page |5[). Such a degeneracy could be 
ascribed for example to a flavor symmetry, the parameter 5 would then characterize 
its breaking. In the simplest case, the flavor group Gf is taken as a Z 2 and the RHNs 
have different parity Z 2 assignments, i.e. Ni ~ odd (even) and N 2 ~ even (odd) 
under Z 2 . The flavor symmetry is broken by the vev of the odd field if). Restricting 
to the 12 block, the resulting mass matrix for the RHNs is 

with 5/2 = (if)) /A. The diagonalization of the mass matrix yields two quasi-degenerate 
RHNs with a mass-splitting 5M . 



3.2 Leptogenesis from soft SUSY breaking A-terms 

In the traditional leptogenesis scenario, the contributions of the soft SUSY breaking 
A-terms to the CP parameter e are usually neglected. Indeed, SUSY breaking will 
induce the following A-terms 

£ so ft = AijmspLiNjHu + h.c. . (14) 

Let us consider the following vertex diagrams, where in the tri-scalar vertex we put 
the A-term contribution from Eqt ( jTjj ) instead of the standard SUSY one. 




k s 



N, 



H 



(a) 



N 



H/ 
- -< 

Lk ^ 



Ni 



H 



(b) 



Figure 1: SUSY breaking ^4-term contributions to the CP parameter e. 

Estimating the contribution of the SUSY soft breaking A-terms to e and com- 
paring it to the standard SUSY one for each of the two considered diagrams, we 
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obtain 

c soft 

£ {a) |/li m 3/2 • x / 1k n 

iSUSY ~ |A|^-sm5 S oft, (15) 

where 5 so f t is an effective soft CP phase. From (pT5|) , we see that in the conventional 
leptogenesis scenario, where the mass of the lightest RHN is Mi ~ 10 10 GeV, e so f t is 
suppressed with respect to e SUSY at least by a factor of 10~ 7 . However, in our scenario, 
where Mj ~ ^3/2, the CP asymmetry parameter e so f t is no more suppressed. It can 
even dominate the over the SUSY contribution depending on the value of the soft 
parameters. This means that CP violation may completely originate from the soft 
SUSY breaking sector, like in the Affleck-Dine case. However, besides enhancing the 
amount of CP violation, the soft SUSY breaking interactions could bring the RHNs 
decay at equilibrium, erasing considerably the produced lepton number. A more 
accurate analysis, requiring the integration of Boltzmann equations, is necessary to 
reach a firm conclusion. 

Finally, it is worth noticing from Eqt. (0) that gravitinos could no more constitute 
a sizable amount of dark matter in our scenario. Indeed, f2 3 / 2 h 2 = 0.01 — 1, requires 
the gravitino to be lighter and/or the gluinos masses to be heavier. 



4 Affleck-Dine leptogenesis with TeV scale RHNs 



Now, we turn to investigate the Affleck-Dine mechanism [JTT], [T^] in the presence of 
TeV scale RHNs. Consider the LH U MSSM flat direction given by 2 

1 f*\ H.. = ±(° 



L -7ll°J- "»=VlUJ (17) 

This flat direction is lifted by the non-renormalizable operator W^nr = A(L H u ) 2 /M = 
A<£> 4 /4M. This operator can be generated via the see-saw mechanism when integrating- 
out the heavy RHNs. The evolution of the scalar condensate <p in the expanding 
background is dictated by the classical equation of motion 

»> + 3 ** + ^ = (18) 

where V(<p) is the full potential, including the soft masses, the Hubble induced masses 
and the A-terms (both from SUSY breaking and the Hubble induced ones)Q. 

,2 



V{<p) = (mi /2 -c H H 2 M z + a H H— 



2 



3 



The factor y2 is necessary to have a canonical kinetic term for ip (The Kahler potential is 



K = H„Hl+LL^ = tpipi 



Here, we are simply ignoring thermal effects 
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4 I \ 1 2 

+a m m 3/2 ^ + h.c. + U_M6. (19) 

The constants ch and an depend on the detailed structure of the Kahler potential. 
In particular the sign of c# is crucial for the validity of the AD scenario. We as- 
sume throughout the paper that it is positive (c# > 0). The evolution of the scalar 
condensate follows three phases. During inflation, when H ^> m 3 / 2 the field (p is 
over-damped and it settles away from the origin at a distance 



\<Po\ * ( JL ]x^) ■ ( 20 ) 



From the last equation, one sees that tp is displaced farther as the neutrino Yukawa 
coupling A is smaller. That is why Lj in Eqt. (O) is usually chosen as the neutrino 
with the smallest Yukawa coupling, L\ say. When H ~ ^13/2, the A-terms enter 
into play and the condensate begins to oscillate. In general, when taking into account 
thermal effects, the condensate begins to oscillate when the decreasing expansion rate 
reaches a certain value denoted H mr , determined when the thermal contributions are 



taken into account [ 35] . At later times when H <C 777.3/2, the lepton number i 



is 



essentially conserved. The evolution of the lepton number, ni defined as 

n L = -((f*(p - (f(f*), (21) 

follows the equation 

ri L + 3Hn L = lm (22) 

The generated lepton asymmetry can be approximated by integrating the equation 
(f22|). This gives 

n L « ^l m (a m ^)t (23) 

In a matter dominated Universe, the expansion rate scales with time as H = 2/31 
Plugging this into the last equation, we get 



s 12 \H osc J \ M, J \MJ |A| 



12 



where M* = Mpi anc k/A/87r ~ 2.4 x 10 18 GeV is the reduced Planck mass and we have 
dropped constants of 0(1). The effective CP- violating parameter S e s is defined as 

5 eff ~ sin (4 arg (p + arg a m ) (25) 

Now, specializing to the low scale see-saw models |22], ^j], where Yukawa couplings 
come-out naturally suppressed as A ~ |Y cfT | 2 ~ m 3 / 2 /M*, we get 



n L 1 / T] 



RH 

12 \ H osc 



S cS (26) 
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Usually the effective CP- violating parameter is assumed to be maximal i.e. £ e ff — 1. 
In our case there is no need to do so, since the reheating temperature can be as 
low as m 3 / 2 so TeV mass RHNs are produced thermally, while the gravitinos are 
not. Typically, the condensate begins to oscillate at H osc > m 3 / 2 . In the extreme 
case when Trh — -f^osc — m 3/2, only a small amount of CP is sufficient to reproduce 
the observed value, namely 5 c g ~ 10~ 9 — 1(U 10 . Up to now, we did not specify the 
transmission mechanism of SUSY breaking. We just assumed that some hidden sector 
will produce the soft breaking scalar masses and A-terms. In the gravity-mediated 
scenario the A-terms are known to be of the form am 3 / 2 W + h.c. This means that 
a m = b m \ and a H = 6#A, where bfj, b m ~ 0(1). In this case, the resulting lepton 
asymmetry is given by 

where now the CP violation parameter is defined as 

5 cS ~ sin (4 arg y? + arg a m + arg A) (28) 

For the typical values H osc ~ m 3 / 2 , T RH ~ 10 9 GeV and 5 cff ~ 0(1), we obtain 
the right value for the lepton asymmetry. 

5 Conclusions 

To conclude, motivated by the potential conflict between the gravitino overproduction 
bound and the high reheat temperature required to produce RHNs thermally, we 
investigated the baryogenesis through leptogenesis scenario in the presence of low 
scale RHNs. We have seen that, in such a scenario the Yukawa couplings of RHNs 
have to be suppressed, in order to give rise to acceptable light neutrino masses. This 
suppression proved to be useful for many purposes, in particular in satisfying the 
out-of-equilibrium condition in the FY scenario. Due to this suppression, however, 
the resulting CP was too small. We used two different mechanisms to enhance the 
CP parameter: the degeneracy of RHNs and soft A-terms. In the latter case, the 
necessary CP violation may come entirely from the soft SUSY A-term. We also 
considered leptogenesis via the LH U flat direction. We have seen that for generic 
SUSY breaking scenarios, AD leptogenesis with low scale RHNs is possible, though 
with reheat temperatures higher than TeV. 
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